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Quantum Monte-Carlo method without negative-sign problem 
for two-dimensional electron systems under strong magnetic fields 
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The quantum Monte-Carlo method is applied to two-dimensional electron systems under 
strong magnetic fields. The negative-sign problem involved by this method can be avoided 
for certain filling factors by modifying interaction parameters from those of the Coulomb in- 
teraction. Our techniques for obtaining sign-problem-free parameters are described in detail. 
Calculated results on static observables are also reported for Landau level filling v = 1/3. 
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It has been revealed that a two-dimensional electron 
system in a strong magnetic field has various ground- 
state phases. In particular, the existence of incompress- 
ible liquids gives rise to the fractional quantum Hall ef- 
fect 1-3 and has attracted much theoretical and experi- 
mental interests to this so-called quantum Hall (QH) sys- 
tem. Rich quantum phases realized in the QH system are 
attributed to Coulomb interactions between electrons, 
because the kinetic energy is quenched by the Landau 
quantization in a strong magnetic field. Thus the ground 
state involves strong electronic correlations and is diffi- 
cult to consider by use of mean field theory. 

So far several many-electron states have been proposed 
theoretically for the ground state in each phase. The 
Laughlin state, 4 Jain state, 5 and Pfaffian state 6 success- 
fully account for incompressible liquid phases in the QH 
system. On the other hand, the composite-fermion liquid 
states 7 ' 8 and anisotropic charge-density- wave states 9 are 
plausible to compressible phases at the half filling of the 
lowest and higher Landau levels, respectively. 

When the properties of such strongly-interacting sys- 
tems are investigated, numerical studies can provide a 
variety of significant information. The exact diagonal- 
ization method of the Hamiltonian matrix has been fre- 
quently used to inspect the various properties of the QH 
system 10 so far, but its applicability is limited to small- 
size systems. Recently the density matrix renormaliza- 
tion group (DMRG) method has been applied to the 
study of the ground-state properties of the QH system. 11 
The system size permitted in the DMRG method is much 
larger than that in the exact diagonalization study, but it 
is difficult to study dynamical quantities by the DMRG 
method. 

In the present study, we propose an application of the 
quantum Monte-Carlo (QMC) method to the QH sys- 
tem. This is because the QMC method can investigate 
not only static but also dynamical quantities in large- 
size systems. This method can be applied for both zero- 
and finite-temperature cases and has ever been used for 
a variety of quantum many-body systems. However, the 



QMC method is accompanied by the negative-sign prob- 
lem, which sometimes makes us away from investigat- 
ing physical properties. Thus we first need to resolve the 
negative-sign problem in order to make the best use of 
the present method. In this letter, we present one solu- 
tion to the negative-sign problem in the QH system. 

This letter is organized as follows. We first express 
physical quantities in the QH system in terms of the 
auxiliary-field path-integral. Then we describe how to 
avoid the negative-sign problem in the present system. 
Finally some results of QMC calculations are presented. 

The system studied here is composed of interacting 
electrons confined on a spherical surface. 12 It is assumed 
that a magnetic monopole is located at the center of 
the sphere. The monopole induces a uniform magnetic 
field on the surface. The number of flux quanta diverg- 
ing from the monopole is denoted by 2s hereafter (2s: 
integer). Single-electron states on the sphere are spec- 
ified by the Landau level index and the z-component 
quantum number, m, of angular momentum (AM). We 
consider the strong magnetic-field limit and restrict the 
single-electron states to the subspace of the lowest Lan- 
dau level (LLL). Then a single-electron state is specified 
by only m which ranges from — s to s, and the number of 
degenerate single-electron states amounts to N s = 2s +1. 
The spin degrees of freedom are neglected for simplicity. 

In the LLL approximation, the kinetic energy is con- 
stant for QH systems containing the fixed number of 
electrons. The remaining term of the Hamiltonian comes 
from Coulomb interactions between electrons. We ex- 
press the interaction Hamiltonian in a bilinear form 
of the density operator to introduce the Hubbard- 
Stratonovich transformation needed for QMC calcula- 
tions. 

The density operator is defined using 
the Clebsch-Gordan coefficient as pkn = 
Y,tn 1 ^ s T,L ? =-s( KN \ sm ^ sm 2) a L 1 a m2 , where a m 
is the annihilation operator of an electron with AM 



quantum number m and a m = (— 1) 



is the time 



reversal of a m . The interaction Hamiltonian is written 
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in terms of the density operator as 



H 



j 2s K 

ozH XK Yl i-^) K+N PKN PK-N + C0P00- (1) 



K=0 N=-K 



Here e = - Ejlot- 1 ) 25 "^ 2 - 7 + l)K//(2\/27TT) and 
the coefficient %x is given in terms of the Haldane 
pseudopotential, 13 Vj, as Xif = ^2j—oTkjVj. The ma- 
trix T is explicitly written using Wigner's 6j symbol as 

Tkj = T,%o(-}) 2 '- J+K+1 (2J + !){'.&}■ 

The interaction Hamiltonian can also be written in 



another form as H = \ EjLo Vj Em= 



-./ A\ m Aj M , 



where Ajm is the pairing operator given by Ajm 
Em 1= - S El 2 =- s ( Jiw l sm i; sm 2 )a mi a mr Since a m sat- 
isfies the anti-commutation relation, Ajm = for even 
2s — J. Thus the terms with Vj for even 2s — J are un- 
physical and have no influence on physical quantities. 
Therefore these unphysical pseudopotentials can be used 
as tunable parameters to improve the efficiency of QMC 
calculations. 

The imaginary-time evolution operator e~@ n is decom- 
posed to the product of imaginary-time slices e~ A/3W (/?: 
inverse temperature, A/3 = /3/N t , N t : Trotter number). 
By introducing the Hubbard-Stratonovich transforma- 
tion for each slice, the evolution operator e~ ,3W can be 
written in terms of a linearized Hamiltonian given by 



Her) 



j 2s r K 

~ X! VKX-K /-I 2 ( a *KNPKN + VKNPKn) 



K=l 



N=l 



+ a K0PK0 + &K0PK0 



+ oXo(Poo) 2 + £oPoo-(2) 



Here pkn — (— ^-) K Pkn-> g kn is an auxiliary field for 
a mode (K, N) introduced by the Hubbard-Stratonovich 
transformation, and t]k is a numerical factor which takes 
1 for non-negative xk and i for negative one. We remark 
that the terms, Xoo(poo) 2 / 2 and eopoo, in Eq.(2) are con- 
stants because poo is proportional to the fixed number of 
electrons. 

The partition functions, Z = Tr[e _ ^ w ] for canoni- 
cal ensemble and Z = (ij)\e~ P n \ij)) for zero tempera- 
ture, are represented in terms of auxiliary-field path- 
integral as Z — J T><t((<t). Here the weight func- 
tion is defined by ((a) = Tr[U((3; 0)] for canonical 
ensemble and by ((a) = (i[)\U((3; 0)|V>) for zero tem- 
perature, U(P; 0) in the weight function is defined by 
U((3; 0) = e - A W- (Nt, )---e- A ' 3 ' l ( CT<1) ), and |V> in the 
zero-temperature formalism is an arbitrary state not or- 
thogonal to the true ground state. The expectation value 
of an observable O is written in the path-integral rep- 
resentation as (O) = JVa(0) a C(a)/Z, where (<D) a = 
Tr[0[/(/?;O)]/C(cr) for canonical ensemble and (0) a = 
(i/)\U(0;0/2)OU(J3/2;O)\iI)}/C(<t) for zero temperature. 

The auxiliary-field path-integral is evaluated by means 
of the Metropolis-Monte-Carlo technique. However, the 
weight function, CC "); can be negative for certain 
configurations of auxiliary fields. Thus let us regard 
P(a) = ICC ") I as a positive-definite probability distri- 
bution. The expectation value is then written as (O) = 
JVa[(0) a P(<T)/fV*Z{v)P(v), where = 



C( (T )/IC( cr )l is the sign of C(ct). If the expectation value of 
£(cr) is vanishingly small, Monte-Carlo evaluations of (O) 
become very unstable. This is the so-called negative-sign 
problem. Thus, in order to make a precise evaluation, it 
is desirable that ((a) is always positive. 

For detailed discussion on the negative-sign prob- 
lem, let us move to the matrix representation of 
the QMC method. 14 In order to exploit the time- 
reversal symmetry, we first define fermionic opera- 
tors, a m , by a m — a m for m > 0, and a. m — 
(— l) s ~' m a m for m < 0. By introducing a matrix M(<r) 
as -A(3h(a) = El^-s^WaL"". we define 
U(/3;0) = e M( CT <N '»)... e M(<r( 1 ))_ The weight function 
in the canonical ensemble is produced by operating 
the particle-number projection to that in the ground- 
canonical ensemble: 15 



((a) = ]T e -2T«v.(VJV.) d et[l + e 2 ™(™ u(/3; 0)], (3) 



i=i 



where N e is the number of electrons. In the zero- 
temperature formalism, the weight function is given by 
C(<r) = dct[V T U(/3; 0)V], where an arbitrary state is 
written in terms of a matrix V and electron vacuum |0) 



as 



W) = (£m=-s °L V ml) ■ ■ ■ (Em=-s "L^ 



.We 

remark that M and U are square matrices of dimension 
N s , while V is a N s x N e rectangular matrix. 

For the zero-temperature formalism, the negative-sign 
problem is overcome under the following conditions: 16 

(i) 2s is odd and N e is even. 

(ii) xk is non-negative for K = 1, 2, • • • , 2s. 

By the condition (i), the dimension, N s , of matrix M 
becomes even. The condition (ii) yields r\K = 1 in 
Eq.(2) and then M satisfies M_ mii _ m2 



^mi,mj & nd 



M_ 



-M* 



-mi,m 2 - _J "mi,-m 2 for positive mi and m 2 . There- 
fore, by arranging AM indices, m, in the matrix repre- 
sentation as m = s, ■ ■ ■ ,1/2, — s, • • • , — 1/2, M has such 
a form as 



M = 



A 

B* 



B 

A* 



(4) 



where A and B are square block matrices of dimension 
N s /2. When matrices Mi and M2 are of the form of 
Eq.(4), the product MiM 2 is also of the form of Eq.(4). 
Thus U(/3;0) also has the form of Eq.(4). When wc 
choose a matrix V that satisfies V- 



the following property is also obtained: 



,N e /2+n 



VZ,n and 



V£ lNe , 2+n for positive m and 1 < n < N e /2, 



V t U(/3;0)V = 



P 



Q 

p* 



(5) 



where P and Q are square block matrices of dimension 
N e /2. Then, if A, is an eigenvalue of VtU(/3; 0)V, its 
complex-conjugate is also shown to be its eigenvalue. 
Thus we obtain ((a) = Ylf^ 2 A,A| > 0. 

We can avoid the negative-sign disaster also in case of 
the canonical ensemble by realizing the following condi- 
tions: (i') 2s is odd, and (ii) xk > for K = 1, • • • , 2s. 
We note that even N e is not required in (f). Since the ma- 
trix M and U(/3; 0) then have the form of Eq.(4) as in the 
zero-temperature formalism, the weight function in the 
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ground-canonical ensemble, det[l + U(/3; 0)], always be- 
comes non-negative. The weight function in the canonical 
ensemble, which is obtained by operating the particle- 
number projection to that of the ground-canonical en- 
semble [see Eq.(3)], is expected to be almost always pos- 
itive. We have numerically confirmed that this is the case. 

The next task is to inquire which systems satisfy the 
negative-sign-free conditions. The condition (i) [(i') for 
the canonical ensemble] restricts the Landau-level fill- 
ing accessible by the QMC method. For instance, the 
v = 1/m Laughlin state (m: odd), characterized by 
2s = m(N e — 1), can be studied by the QMC method 
if N e is even. The condition (ii) brings about strong con- 
straints for electron-electron interactions accessible by 
the QMC method. In fact, this condition is not satisfied 
by the Coulomb interaction within the LLL. 

As mentioned above, the Haldane pseudopotentials, 
Vj, are classified into the physical components (2s — 
J: odd) or unphysical ones (2s — J: even). We ob- 
tain coupling constants xk favorable to QMC calcu- 
lations by using the unphysical ones as tunable pa- 
rameters, with minimizing the difference of the phys- 
ical ones between our model interaction and the 
Coulomb interaction. Then the model pseudopoten- 
tials are given in terms of obtained xk by Vj = 
YJk=o t JkXk, because J2k t JkT K j> = 5jj>. We define 
the quantity to be minimized by /(— xa, Xii ' • • > X2s) = 

E 2s -J:odd K J (J2k=o T JKXk - where Vp is the 

pseudopotential of the Coulomb interaction and kj con- 
trols the priority of each pseudopotential component. We 
note that the arguments of / should be non-negative, 
which agrees with the condition (ii). Although the first 
argument, — Xoi is n °t referred by the condition (ii), it 
must be positive in order to make a repulsive pseudopo- 
tential Vq = J2k=o TokXK > under the conditions that 
Xk > for 1 < K < 2s and T 0K < for all K. Such 
a set of xk can be obtained by the linear programming 
method, which minimizes / over its non-negative argu- 
ments imposing the inequalities Vj > Vp (Vj < Vp) for 
non- negative (non-positive) kj. 

A simple condition for the linear programming method 
is given by kj = 1 for all of odd 2s — J. Then physi- 
cal pseudopotentials are determined with equal weights 
and become larger than those of the Coulomb interac- 
tion, respectively. The coupling constants xk obtained 
in this way vanish except for the cases of K = and 
K = 1. The Hamiltonian for this interaction is written 
as H = -±xoN2/N s +6 X i£ 2 /[(N s -l)N s (N s +l)], where 
C denotes the total AM operator. The eigen-energies of 
this Hamiltonian depends only on the total AM quan- 
tum number. Therefore the ground state is energetically 
degenerate and the excitation spectrum is gapless. 

Figure 1 shows the pseudopotentials of the Coulomb 
interaction and two negative-sign-free interactions. For 
the interaction I whose pseudopotentials are shown by 
cross symbols, we set the priority parameters as kq = — 1 
and kj — for J ^ 0. Namely we minimized the differ- 
ence, | Vb — Vp\, between the longest-range components to 
mimic the long-range nature of the Coulomb interaction. 
As the inequalities imposed in the linear programming 



^0-5 
-3 0.4 
0.3 



« i 


i 

i 


i 

- ) 

< 


< » 


i L 

i ; 


: 


: 




< 






— Coulomb 
X model 1 
A model 11 

" 1 1 " 1 1 " 1 1 " 




< 











O(N^^00O(N^r^C0O 



Fig. 1. Pseudopotentials for negative-sign- free interactions and 
the Coulomb interaction within the lowest Landau level with 
N s = 2s + 1 = 22. Only physical components for odd 2s — J are 
shown. The pseudopotentials of negative-sign-free model inter- 
action named I are obtained by the linear programming method 
and are shown by cross symbols. The pseudopotentials of in- 
teraction II shown by closed triangles are obtained by applying 
uniform shifts and rescaling to those of interaction I. 

method, we used the following ones: V2S-1 > Vp s _ l (the 
shortest-range component) and Vj < Vp for J ^ 2s — 1. 
This is because the Laughlin state favors prominently 
large short-range components. Pseudopotentials of model 
interaction I agree well with those of the Coulomb inter- 
action on short-range components, but do not on long- 
range ones. In fact, the long-range components shift 
downward from those of the Coulomb interaction and 
eventually become even negative. 

In order to make improvements on the interaction I, we 
perform the following operations on its pseudopotentials: 
(i) uniform shift of pseudopotentials, Vj — > Vj + \Vq\, (ii) 
energy rescale, Vj -> Vj x (Vg^ - V?)/(Vl,_ x -Vj), 
and (iii) another uniform shift of pseudopotentials, Vj — > 
Vj + Vp. As a result of these sequential operations, 
we obtain another interaction named II, whose pseu- 
dopotentials are given by V} 1 — Vp + (Vj - Vq) x 
(Vp s _ 1 - V r c )/(V A 2 I S _ 1 - Vj). We note here that uniform 
shifts of pseudopotentials make no influence on the elec- 
tronic states, because the change, Vj — > Vj + v for all 
J, involves just a constant term in the Hamiltonian as 
^vJ2 JM A JM Aj M — ^vN e (N e — 1). Furthermore these 
operations do not break the negative-sign-free character 
of xk i because uniform shifts of Vj bring about only the 
change in xk=o and energy rescale keeps the sign of xk- 
Thus the interaction II is also free of negative-sign. 

The pseudopotentials of interaction II are shown by 
closed triangles in Fig.l. The energy scales of interac- 
tion II and Coulomb interaction are the same, because 
V2S-1 — Vb gives the energy scale. The character of in- 
teraction II is similar to that of the Coulomb interac- 
tion on short- and long-range components, while the dif- 
ference between pseudopotentials of the interaction II 
and Coulomb interaction is not so small on middle-range 
ones. In spite of this disagreement on middle-range ones, 
it is shown later that the interaction II leads to the 
ground-state properties qualitatively similar to those of 
the Coulomb-interacting system. 

Figure 2 shows the results of zero-temperature QMC 
calculations for the QH system with 2s = 21 and N e = 8 
[y = 1/3). These calculations were done by using the 
negative-sign- free interaction II in Fig.l. Calculated re- 
sults for three imaginary-time slice widths, A/3 = 1/4, 
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Fig. 2. Expectation values of (a) energy per particle and (b) total 
angular momentum by the zero-temperature QMC calculations 
of the QH system with 2s = 21 and N e = 8 {v = 1/3). The 
unit of energy is e 2 /lg (lg: the magnetic length), j3 the inverse 
temperature, and A/3 is the width of imaginary-time slice. The 
energy contribution by the positive-charge background is also 
taken into account. 



1/8, and 1/16 [1/ (e 2 /Ib)} are almost the same. Hence 
the numerical errors introduced by the Suzuki- Trotter 
approximation are not so important for A/3 < 1/4. For 
the inverse temperature (3 > 60/(e 2 /Zs), both expecta- 
tion values of energy and angular momentum are almost 
saturated. The energy gap for the interaction II is found 
to be AE ~ 0.09(e 2 /^s) by the exact diagonalization 
method. The inverse temperature needed for convergence 
is so high (P > 6/AE). This seems to be due to a high 
density of states for E > AE. 
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Fig. 3. Expectation values of (a) energy per particle and (b) 
ODLRO correlation function Cs(r). Results by the exact diag- 
onalization and finite-temperature QMC calculations are shown 
for the QH system with 2s = 21 and N e = 8 {y = 1/3). r arc 
denotes the arc-distance between r and the north-pole on the 
Haldane spherical surface. 



Figure 3(a) shows the /3-dependence of energy expec- 
tation value per particle in the canonical ensemble. The 
results by the exact diagonalization and QMC calcula- 
tions are in good agreement. For > 60/(e 2 /Ib), the 
expectation values of energy per particle are almost satu- 
rated as well as in Fig. 2 (a). Since the existence of the off- 
diagonal long-range order (ODLRO) is one of the most 
remarkable properties of the Laughlin state, 17 we calcu- 
lated the ODLRO correlation function in the v = 1/3 
QH system with interaction II. The correlation function 
is defined 18 by C 3 (r' - r) = (A t (r')^4(r)), where ^4(r) is 



the annihilation operator of a composite boson (an elec- 
tron with three flux quanta attached to it). Figure 3(b) 
shows C 3 (r) for (3 = 8 and 64 [l/{e 2 /l B )} obtained by 
QMC calculations as well as that at T = obtained by 
the exact diagonalization. At zero temperature, Cq(t) re- 
mains finite for large separation r. This result suggests 
that the ground state for the interaction II should have 
the ODLRO. For finite temperatures (3 < 8, the ODLRO 
correlation rapidly vanishes with increasing separation 
r, while it survives at j3 = 64 ~ 6/AE as in case of 
T = 0. This is because the contributions of excited states 
to Cs(r) is negligible for temperatures quite lower than 
AE. 

We comment on the validity of negative-sign-free in- 
teraction II as a substitute for the Coulomb interaction. 
By the exact diagonalization method, we confirmed that 
the overlap between the two ground states for the inter- 
action II and Coulomb interaction is not large (~ 0.35) 
for v = 1/3. However, the existence of the ODLRO in 
the ground state and that of energy gap above it sug- 
gest that an imperfect Bose condensation of composite 
bosons exists for interaction II, although it is not so per- 
fect as that in the Laughlin state. Thus we believe that 
the ground state for interaction II contains the essential 
properties of the Laughlin state. 
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